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ABSTRACT 

In recent years, we have come to recognize the widespread importance of large-scale winds in the 
lifecycle of galaxies. The onset and evolution of a galactic wind is a highly complex process which 
must be understood if we are to understand how energy and metals are recycled throughout the 
galaxy and beyond. Here we present 3D spectroscopic observations of a sample of 10 nearby galaxies 
with the AAOmega-SPIRAL integral field spectrograph on the 3.9m AAT, the largest survey of its 
kind to date. The double-beam spectrograph provides spatial maps in a range of spectral diagnostics: 
[Oiii]5007, H/3, Mg6, Nao, [Oi]6300, Ha, [Nii]6583, [Sii]6717, 6731. We demonstrate that these 
flows can often separate into highly ordered structures through the use of ionisation diagnostics and 
kinematics. All of the objects in our survey show extensive wind-driven filamentation along the minor 
axis, in addition to large-scale disk rotation. Our sample can be divided into either starburst galaxies 
or active galactic nuclei (AGN), although some objects appear to be a combination of these. The total 
ionizing photon budget available to both classes of galaxies is sufficient to ionise all of the wind-blown 
filamentation out to large radius. We find however that while AGN photoionisation always dominates 
in the wind filaments, this is not the case in starburst galaxies where shock ionisation dominates. 
This clearly indicates that after the onset of star formation, there is a substantial delay (> 10 Myr) 
before a starburst wind develops. We show why this behavior is expected by deriving "ionisation" and 
dynamical timescales for both AGNs and starbursts. We establish a sequence of events that lead to 
the onset of a galactic wind. The clear signature provided by the ionisation timescale is arguably the 
strongest evidence yet that the starburst phenomenon is an impulsive event. A well-defined ionisation 
timescale is not expected in galaxies with a protracted history of circumnuclear star formation. Our 
3D data provide important templates for comparisons with high redshift galaxies. 

Subject headings: Galaxies: individual (NGC 253, NGC 1365, NGC 1482, NGC 1808, NGC 3628, 
NGC 5128, Circinus, NGC 6240, NGC 6810, IC 5063) 



1. INTRODUCTION 

Recent observations reveal that outflows are important 
in the lifecycle of galaxies. First, Tremonti etal. (20041 
find evidence for chemical enrichment trends throughout 
star-forming galaxies over three orders of magnitude in 
stellar mass. Surprisingly, the enrichment trends are ob- 
served to continue all the way down to log(iboi/i0) ~ 4 
( Kirby e^a/.|2008 ) where Lboi is the bolometric luminos- 
ity of the galaxy. Indeed, there is evidence that galaxies 
that exceed the mass of the Milky Way manage to retain 
a large fraction of their metals, in contrast to lower mass 
galaxies where metal loss appears to be anticorrelated 
with baryonic mass. Recent theoretical work has shown 
that a low-metal content in dwarf gala xies is attributabl e 
to efficient metal-enriched outflows ( [Dalcanton 2007) 
Secondly, there is now strong evidence tor large-scale ou 
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consistent with the action of winds (Cen &: Ostriker|1999 



Madau e^aL||2001| |Ryan- Weber etai 
2008| . A high proportion of Lyman 



Dave et al. 
reak galaxies at 
4 show kinematic signatures of winds (Erb etal^ 
12006') as do many nearby dwarf starbursts ( Schwar tz H 
Martin 2 0d4| and ultraluminous infrared galaxie s (Mar 
^ i||2006[ |SS& 
N-bod' 



^1 



iwartz, Martin, Chandar aZ.||2006| . 



-body simulations of galaxy formation within cold 
dark matter (CD M) cosmology fail to p r oduce realis 



tic galactic disks (Navarro & White 1994 Steinmetz & 
Miiller 1995| ) . Over-cooimg in CDM-f-hydrodynamic sim- 



ulations, due primarily to a lack of resolution, results in 
too much gas arriving in the centre of the galaxies. But 
a useful aspect of these incomplete models has been to 
highlight the possible role of feedback in shaping galax- 
ies. Vigorous feedback in the early phase of galaxy for- 
mation is a partial solution to the angular momentum 
problem by preventing baryons from losing too much of 
their specific angular momentum through the interac- 
tion with dark matter (Fall 2002 ). Indeed, some authors 
have claimed more realistic disks retaining much more of 
their angular mome ntum when one includ es the action 
of AGN-drive n jets (Robertson e^aL || 2006 ) or starburst 
driven winds ( Efstatlirou||2000[ ) '. 

Our new survey is preernpted by a vast literature across 
a broad spectrum that can be traced back over two 
decades — comprehensive reviews are given by Heck- 
man efaL (1989) and YeiWeux etal. (2005). Much of 
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the early optical work on galactic outflows was lim- 
ited to narrowband images and long-slit spectroscopy 
(e.g. M cCarthy g/. (1987); [Heckman, Armus fc Miley 



(1990); Phillips (1993); Lehnert & Heckman (1996)) al- 



though Fabry-Perot observations have been carried out 
in a few emission lines (e.g. Bland & TuUy (1988); Ce- 
cil (1990)). The quality of UV, x-ray, infrared and ra- 
dio continuum observations has greatly improved in re- 
cent years (e.g. Hoopes et al. 2005; Str ickland et al. 2004; 

Even without 



( |2009| ). 



Grimes, Heckman, Aloisi et al. 

the kinematic separation of compound emission struc- 
tures that is possible via spectroscopy, there are clear 
associations between these bands and the optical maps 
(e.g. Veilleux et 1994; Cecil etaZ. 2002). Historically, 
most winds had been confirmed kinematically from evi- 
dence of line-splitting along the galaxy minor axis, but 
at least one wi nd has been di scovered through its ionisa- 
tion signature ( Veilleux fc Rupke,,2002) . In recent times 
measurement oi blueshitted absorption lines have proved 
powerful in this regard. 

In this new study, we perform wide field integral field 
spectroscopy of a sample of nearby galactic winds. The 
advantages offered by a wide wavelength base line and 
moderate spectral resolutions allow the determination of 
physical parameters for wind outfiow mechanics based 
on ionisation and kinematic signatures e.g. rotating disk, 
high-latitude warm ionizing medium, wind filaments, en- 
trained gas. One well-studied object in particular con- 
vinced us of the potential for wide-field mapping with 
an integral field spectrograph. The projected distribu- 
tion of optical line emission in M82 is exceedingly com- 
plex but with the aid of Fabry-Perot observations. Shop- 
bell & Bland-Hawthorn (1998) were able to kinematically 
separate the front side from the back side of the wind 
flow, and from the diffuse halo of broadened line emission 
see also the recent IF S studies o f Wcstmoqucttc, Smith, 

(i2r~ 

^he 



Gallagher et al. mm ] 



^ Westmoquette, Gallagher, Smith| 

e?a/. (2009])). When this is done, both the ionisation and 
Kinematics exhibit well ordered behavior. The lesson is 
clear: with complete observations, it is possible to sep- 
arate distinct dynamical components. Only then does 
it become possible to determine physical parameters for 
the outflow mechanism. 

In this flrst paper, we present the new observations 
and describe the reduction procedures specific to inte- 
gral field spectroscopy (§[2]). We present emission-line 
flux and line ratio maps ror a wide range of diagnos- 
tics (e.g. ionisation and density). For each galaxy, we 
use the ionisation diagnostic diagrams (IDDs) to delin- 
eate the different sources of ionisation as a function of 
location and ionisation source (see §p| and classify the 
sample targets acco rding to their individual IDDs (§|3.1| 
- starbursts 
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AGN). These diagnostic diagrams 
reveal well-ordered behavior across each galaxy. In §|4] 
we consider the source of ionisation from a theoretical 
perspective and use these considerations to construct an 
evolutionary sequence for galactic winds in section §[5j 
Our conclusions are presented §|6] 

Through this work, we assume a flat concordance cos- 
mology such that Ho=71 kms~^ Mpc^^, riM=0.27 and 
1}a=0.73. 

2. OBSERVATIONS AND DATA REDUCTION 



2.1. Observing set up and procedure 

The AAOmega integral-field spectrograph at the AAT 
3.9m is ideally suited to our study (see Appendix). 
The SPIRAL integral field unit (IFU) is a 32x16 el- 
ement rectangular microlens array coupled via an op- 
tical fi bre feed to the dua l beam A AOmega spectro- 
graph ( Saunders et al. 2004 ; Sharp et al. 2006 ) . Mounted 
at the Cassegrain auxilliary focus, the 0.7" IFU pixel 
scale is well matched to the typical seeing at the AAT, 
and gives excellent sensitivity to extended low surface 
brightness features. The single pointing field-of-view is 
22.4" X 11.2". High accuracy, guide-probe offset guiding 
allows simple and reliable mosaicing of large areas. The 
default position angle of the IFU places the long axis 
in an east-west direction. Observations of a number of 
targets were undertaken at alternate PAs to accommo- 
date the demands of other observing programs under- 
taken alongside the observations reported here. 

Observations were undertaken over a number of ob- 
serving runs as detailed in Table [T] The dual beam 
AAOmega spectrograph provides excellent coverage of 
the spectral features of interest, while maintaining a 
high spectral resolution. We observe with the 1500V 
and lOOOR Volume Phase Holographic (VPH) gratings 
in the blue and red arms respectively. We use a 5700A 
dichroic beam splitter and center the wavelength cover- 
age and VPH blazes at 5000A and 68OOA, resulting in 
2.4 pixel mid-range resolution elements of R~6300 and 
R^5600 (velocity resolutions of sa 35 and ^50 km s~^ 
FWHM). On accounting for the sample redshift range, 
this targets the emission lines H/3 / [O ill] in the blue and 
[O i]/Hq:/[N ii]/[S 11] in the red. Unfortunately the im- 
portant [0 11] doublet at 3727A falls outside of the ob- 
servable spectral range with the chosen settings. We con- 
sider this to be a crucial line for future studies of this 
kind. 

Quartz-Halogen flatfield exposures and CuAr-|-FeAr 
arc lamp exposures are taken at intervals throughout the 
observation to allow fibre tracing on the CCDs (both 
arms of AAOmega use E2V 2kx4k CCDs), fiat field- 
ing and wavelength calibration. Stationed in the west 
Coude room, AAOmega is kept thermally and gravita- 
tionally stable throughout each night's observations in 
order to minimize the overhead on system calibrations. 
The strong night sky oxygen line at 5 5 77 A falls outside 
the observed range and so cannot be used to calibrate 
relative fibre transmission as is the common practice in 
low resolution multi-object spectroscopy. Therefore twi- 
light fiatfield frames are observed in order to account for 
relative fibre-to-fibre transmission variations. 

An assortment of spectrophotometric standard starfj^ 
were observed throughout each observing block, in or- 
der to prepare a fiux calibration solution. The scatter in 
the calibration function derived from the standards dur- 
ing the observations indicates a 10% uncertainty in the 
absolute fiux calibration level. The relative calibration 
with wavelength, a quantity more relevant to line ratio 
determinations, is consistent at the ~2-3% level, as de- 
termined from the comparison of calibration observations 

* Spectrophotometric standard stars were chosen from 

the ESQ spectrophotometric standard ^lis t available at 

|http://www.eso.org/sci/observing/tools/standards/spectra/| 
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TABLE 1 

a summary of aat observations with the 
AAOmega-SPIRAL spectrograph. 



Observing Dates Objects Observed 

May 2007 NGC3628, NGC6810, NGC 5128(CenA) 

October 2007 NGC 1365, NGC 1482, NGC 1808 

July 2008 Circinus, NGC 253, IC 6063 

August 2008 NGC 6810, NGC 1705, NGC 6240 



over the course of tlie survey. 

Each target was observed as part of a multi-pointing 
mosaic. An overlap of 1 or 2 IFU elements (0.7" or 1.4") 
was allowed to enable relative intensity scaling and accu- 
rate mosaicing of the final data sets. Individual science 
exposures were typically 1200 sec. Mosaic sizes vary be- 
tween galaxies, but are typically of the order of six to 
twelve IFU positions, yielding mosaics of the order 20" 
to 60" on a side. 

Since our targets completely fill each IFU pointing, it is 
not possible to generate a sky subtraction solution from 
fibres within each IFU pointing. To achieve accurate 
sky subtraction, a small number of dedicated offset sky 
frames was taken each night. Data was processed in real 
time at the telescope using the 2df dr software package. 
This allowed monitoring of the observing conditions such 
that the mosaic pattern could be repeated or adapted if 
required. The final data analysis was carried out using a 
custom suite of idl routines. 

There are four isolated "dead elements" within the 
SPIRAL IFU. During the reduction process, these are 
replaced with the average spectrum of the four adjacent 
IFU elements prior to mosaicing the data. Individual 
frames were then aligned and mosaiced using telescope 
offset information and a suite of custom written software 
tools. Frames are scaled, based on an iterative compari- 
son of overlap regions in the mosaic, in order to account 
for minor variations in transparency, seeing and exposure 
time. 

2.2. The sample 

We targetted a sample of ten galaxies with known 
galactic winds ( Veilleux ei a/. 2003; Heckmaneia/. 1990; 
Lehnert & Heckman 1996). The details of the sources 
are presented in Table [2] Given that the aim of our 
study is to compare the properties of AGN vs. starburst 
winds, we did not include a "control" sample of non- 
wind galaxies. Integral field observations are extremely 
time consuming and indeed IFU surveys of galaxies have 
rarely been attempted before. Notable examples are the 
SAURON survey of early type galaxies (e.g. de Zeeuw 
etal. 2002) and the SINS survey at higher redshift ( Gen 
zel et aL|[2008l ) 



We limit our survey to southern hemisphere sources 
with bolometric luminosities in the range Lboi — 1 — 10 x 
10^" Lq to ensure the impact of photoionisation by the 
central source (whether starburst or AGN) is detectable 
out to large galactocentric distances. (The physical mo- 
tivation for this statement is given in §4.1.) With a view 
to going after a lower luminosity sample, we attempted 
a preliminary study of the d warf galaxy NGC 1705 ( An- 
nibali, Greggio, Tosi e^a/.][2003; (Meurer, Freeman, Do- 



pita & Cacciari 1992 Hunter, Hawley & Gallagher 1993 1 . 
Our initial observations (not presented here) served to 
emphasize that low mass galaxies typically have lower 
metallicities which render many of the ionisation diag- 
nostics intrinsically weak. A survey sample of a substan- 
tial number of low-luminosity wind sources (^10) will 
require a suitable IFU spectrograph on an 8m class tele- 
scope. 

For each of the ten wind galaxies, the SPIRAL IFU was 
used in a simple tiling pattern to cover most of the fila- 
mentary emission above a surface brightness of 1 x 10""'^^ 
erg cm~^ s~^ arcsec"^ (cgs). The SPIRAL mosaic 
patterns (or footprints) are overlayed on _ffST/WFPC2 
imaging data in Fig. [Tj WFPC2 images have been used 
for all sources with the exception of NGC 1482 for which 
UKST data from the Digitized Sky Survey has been used. 
The HST /WFPC2 data were obtained from the Multi- 
mission Archive at the Space Telescope Science Institute 
(MAST). Fig.[T]emphasizes the need for much wider field 
IFUs for the study of nearby galaxies. While these are 
under development in Europe and Australia, we stress 
the simultaneous requirements of spectral coverage and 
resolution (see Appendix) which are often sacrificed in 
order to achieve a larger spatial format. In this respect, 
we consider the AAOmega-SPIRAL spectrograph to be 
ideal for the proposed science, even with the limited field 
of view. 

2.3. Emission line maps 

Emission line maps are constructed via profile fit- 
ting to each individual spectrum from the data cubes. 
Lines are approximated with single Gaussian profiles in 
the work that follows, and each line is fitted indepen- 
dently, rather than simultaneously, which is possible at 
the SPIRAL resolution in most instances. This simpli- 
fication overlooks effects such as line splitting observed 
in emission from the front and back surfaces of an out- 
flow cone and asymmetric profiles (blue/red shoulders to 
emission lines) from marginally resolved kinematic com- 
ponents. Ultimately a full deconvolution of each data 
cube, informed b y kinematic models for each source, will 
be required (e.g. Westmoquette, Smith, Gallagher et al. 
(j2009|). 

Currently no account is being taken for the underly- 
ing stellar absorption although II/3 absorption is evident 
close to the disk in some cases. The limited wavelength 
coverage, which excludes Balmer lines higher than H/3, 
hampers a proper treatment of the underlying stellar con- 
tinuum. The dominant effect of neglecting stellar absorp- 
tion will be to underestimate the flux in the II/3 emission 
line, suggesting that [O iii]/H/3 ratios may be moderately 
overestimated in regions with a significant contribution 
from stellar continuum. In our present study, Balmcr 
absorption has negligible impact since we are concerned 
with filaments that extend beyond the projected stellar 
disk. 

2.4. Ionisation diagnostic diagrams 

For each object, we construct "i onisation diagnostic 
diagrams" (IDD s) as prescribed by Baldwin , Phillips & 



Terlevichl ( [T98l| ) and (Veilleux fc Osterbrock| ( |1987| 
different families of optical emission lines. An exan 



tor 

optical emission lines. An example, 
constructed from the observations of NGC 1482 (Figs. |2] 
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TABLE 2 

The target sample of southern or equatorial galactic wind sources. 



Object 


RA/Dec (J2000) 


Ptedshift 


Hubble type"'^ 


Spec class^ 


i 


Mb 


LboI^ 


R/Mpc 


kpc /arcsec 


NGC 253 


00 47.6 -25 18 


0.00081 


Sc(X) 


Hii 


86 


-20.02 


2.8 


2 


0.017 


NGC 1365 


03 33.7 -36 08 


0.00546 


Sb(B) 


Hll, Sy2 


63 


-21.26 


9.3 


16.9 


0.100 


NGC 1482 


03 54.7 -20 30 


0.00639 


Sa(P) 


Hii 


58 


-18.89 


1.1 


19.6 


0.121 


NGC 1808 


05 07.7 -37 31 


0.00332 


SO/a(X) 


Hii, Sy2? 


50 


-19.52 


2.2 


10.8 


0.059 


NGC 3628 


11 20.3 +13 37 


0.00281 


Sb(P) 


Hii, liner 


87 


-19.96 


2.8 


7.7 


0.041 


NGC 5128 


13 25.3 -43 01 


0.00183 


SO, Lenticular 


Sy2 


43 


-20.97 


8.5 


4.9 


0.053 


Circinus 


14 13.2 -65 20 


0.00145 


Sb(A) 


Sy2 


65 


-21.23 


9.0 


4.2 


0.014 


NGC 6240 


16 53.0 +02 24 


0.02445 


10 


LINER, Sy2 




-21.30 


8.2 




0.5 


NGC 6810 


19 43.6 -58 40 


0.00678 


Sa(A) 


Hii, Sy2? 


82 


-20.61 


5.5 


25.3 


0.135 


IC 5063 


20 52.0 -57 04 


0.01135 


Sa 


Sy2 




-20.34 


4.2 




0.22 



Note 1: Vizier database a t|http://vizier.u-strasbg.fr/| 

Note 2: NED database at |littp: / /nedwww.ipac.caltech.edu/| a question mark indicates that tlie Seyfert classification has been questioned 
in later studies. 

Note 3: Bolometric luminosities, LboI, in units of 10^" L©, are computed from the modeling of |Buzzoni| | |2005| , with uncertainties of 
typically 10% for different assumed ages of the stellar population over the full range 1-14 Gyr. 



TABLE 3 

Classification of the ionisation mechanism for the 
galactic wind emission. 



Object 


[O III] & 


I [N II] [O III] i 


z [S II] [O III] & [O i] 


Starbursts 








NGC 253 


10:90 


10:90 


10:90 


NGC 1482 


10:90 


10:90 


10:90 


NGC 1808 


20:80 


0:0 


0:0 


NGC 3628 


0:01 






NGC 6810 


0:0^ 


0:0 




AGNs 








NGC 1365 


90:10 


90:10 


100:0 


NGC 5128 


50:50 


0:100 


60:40 


Circinus 


100:0 


60:40 


100:0 


NGC 6240 








IC 5063 


100:0 


90:10 


100:0 



The primary ionisation source in the gala ctic wind-driven fila- 
ments is determined from the IDDs in Figs. |4|7|10|12|15|18|21|24| 
& |28[ For each of the IDDs, we show the percentage ot points 
(Agin :Shoc k) outsid e of the maximal star formation boundary 
p<ewley ef ai.||2001[ l consistent with the fiducial models derived 
trom NGG 1365 and NGC 1482. Reliable fractions are hard to de- 
rive so the numbers shown are indicative. A ratio of 0:0 indicates 
that all emission is contained within the star formation boundary. 
Blank entries indicate that the IDD could not be constructed in 
the given line ratio. 

Note 1: Figs. [12] and |14| clearly delineate enhanced [Nll]/Ho co- 
spatial with line profile splitting in the outflowing gas, consistent 
with the shock signatures seen in NGC 1482 (see §3.1). 



&|3]), is presented in Fig. [4} We hav e included the extreme 
starburst model of Kewley etal. (|2001|) as a common 



starburst model of 
fiducial indicator, 
presented in Figs 



Kewley etal. (j2001| as a common 
Jhe IDDs lor the galaxy sample are 
0[Tg[T§[T8l[2l}|4]&[28l 



We have added two further reierenceTines to each di- 
agnostic plot, which are derived from two objects in our 
sample. The two diagonal lines across the IDDs indicate 
the different tracks traced by the bulk of the data points 
in NGC 1365 and NGC 1482. These two galaxies are the 
clear est examples within our sample of AGN p hotoionisa- 
tion ( Veron, Lindblad, Zuiderwi .ik e^a/. | 1980 ) and shock 
ionisation ( Veilieux & Rupke|2UU2[ ) respectively. As such, 
these tracks constitute our own internal calibration of 
two very different sources of gaseous ionisation and al- 



low for a very precise differential analysis for our galaxy 
sample. This clear distinction arising from the IDDs for 
each galaxy allows us to cleanly separate most of the 
galaxies into one or other category in terms of the dom- 
inant source of ionisation in the outer filaments. 

We emphasize that shock diagnostics can overlap a re- 
gion of the IDDs occupied by so-called LINER objects, 
some of which may be excited by dilute photoionizing 
sources rather than by shocks. The origin of LINER ion- 
ization is still somewhat controversial. For all of our ob- 
jects, the strength of the [Nil], [SII] and [OIII] hues, and 
sometimes [01] lines, argues against "LMER" ionization 
so we do not consider it here. 

3. THE GALAXY SAIVIPLE 

We provide an overview of each galaxy based on an 
extensive literature review with particular emphasis on 
recent work. Our initial classification below in terms of 
starburst or AGN is taken from Table 2. We start with 
the canonical objects NGC 1482 and NGC 1365 in each 
category respectively. 

3.1. Starbursts 
NGC 



Hameed & Devereux (1999) first reported "filaments 
and chimneys ot ionised gas" in NGC 1482. IMore 
recently, tunable f i lter o bservations of NGC 1482 led 
Veilieux & Rupke ( 2002 ) to consider it the prime ex- 
ample ot a "shock excited, limb-brightened conic wind 
structure" . No compelling evidence for an AGN has yet 
been reported in the literature at any wavelengt h and 
thus NGC 1482 is classified as a typical Hii galaxy ( Kew- 
ley e^aLpOlj) . 

Our SPIRAL IFU observations of NGC 1482 are pre- 
sented in Figs. |2l-|4l The conic outflow structure seen in 



the emission-line gas ( Veilieux fc Ru pke 2002) and dif- 
fuse x-ray emission ( Strickland et al. 2004a) is evident 
in the SPIRAL integrated line flux maps (Fig. [2|. The 
enhanced line emission in the forward facing component 
of the limb brightened bi-cone is clearly visible above 
the disk (particularly towards the bottom of the [Nil] 
and [Oiii]maps of Fig. [2]) with the counter cone visible, 
subject to dust attenuation, below the inclined disk in 
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TABLE 4 

Spatial extent and brightness of the emission-line gas. 



Object 




Extent 


Deprojected 


/x(Ha) @ 1 kpc 


L(Ho) @ 1 kpc 




(10^10 Lo) 


(arcsec) 


(kpc) 


(ergs s^^ cm~^ arcsec"^) 


(ergs s^^ arcsec"-^) 


NGC 253 


2.8 


30.5 


0.5 


3.8e-16 ± 4.6e-16 


5.4e+35 ± 6.5e+35 


NGC 1365 


9.3 


16.4 


1.8 


2.6e-16 ± 1.2e-16 


1.46+37 ± 6.5e-|-36 


NGC 1482 


1.1 


15.5 


2.2 


8.9e-16 ± 4.4e-16 


6.6e+37 ± 3.3e+37 


NGC 1808 


2.2 


22.5 


1.7 


1.6e-16 ± 5.1e-17 


3.7e+36 ± 1.2e+36 


NGC 3628 


2.8 


28.0 


0.4 


5.2e-17 ± 2.2e-17 


1.7e+36 ± 6.9e+35 


NGC 5128 


8.5 


36.0 


2.8 


4.7e-16 ± 2.2e-16 


6.8e+36 ± 3.2e+36 


Circinus 


9.0 


24.6 


0.4 


1.9e-16 ± 8.5e-17 


1.5e+36 ± 6.6e+35 


NGC 6810 


5.5 


12.0 


16.4 


5.5e-16 ± 2.4e-16 


4.7e+37 ± 2.1e+37 


IC 5063 


4.2 


6.0 


1.4 


1.3e-15 ± 5.3e-16 


3.1C+38 ± 1.3e+38 



The spatial extent of emission-line gas is estimated from measurements of the Ho and [N ll] emission along the minor axis of each galaxy. 
The point at which line fitting becomes unreliable is determined for each mosaic. For the default observing mode, this corresponds to a 
line surface brightness of 2 X 10~^^ cgs. The physical extent within the mosaics is corrected for the galactic inclination given in Table |2| 
In the final two columns, we give the mean Ha surface brightness and the mean Ha luminosity at 1 kpc (with la uncertainties), or at the 
physical limit of the gas if this is less than 1 kpc. These values are not corrected for internal dust extinction. 



most emission lines. Line profiles fi'om the near and far 
side of the extended filaments show clear evidence of line 
splitting at the level of ^130 km s~^. 

Strong line emission is seen throughout the galactic 
disk with the base of the outfiow cone coincident with 
the outer boundary of strong star formation (as mea- 
sured by the gradient in the Ha emission along the disk) , 
indicative of a wind launched above (and below) the 
plane of the disk from an extended starburst region. A 
sharp gradient in the electron density, as measured by 
the [S ii] : [S ii] ratio (Fig. Isl) is also coincident with the 
boundaries of the base of tne outflow region, suggestive 
of a clearing of the inner disk by the action of the star- 
burst wind. 

Under an alternate intensity scaling, the Ha integrated 
flux map does show a compact nuclear region at the 
limit of the spatial resolution of the SPIRAL data. The 
nucleus appears to be res olved at mid-IR wavelengths 
( Siebenmorgen et aL|2008" ) and with no spectroscopic ev- 
idence ot an AGiN the emission is likely that of a compact 
star forming region. 

The line ratio diagnostics for NGC 1482 are shown in 
Fig.llj The limited sensitivity to [Om] emission, partic- 
ularly in the counter flow region to the north of the disk 
(presumably hampered by significant dust attenuation) 
restricts the region that can be analysed in this man- 
ner. As for all of our sample, most of the data points are 
bunched to the right of the IDD, consistent with metal 
enriched gas in the nuclear regions of massive galaxies. 
The body of the disk emission i s consistent w i th hig h 
metallicity star formation (Z>1.0, Kewley etal. (2001)). 

The emission from the bi-cone ot INGG 1482 is com- 



etal. (2004b). 



patib le with shock- induced ionisation ( Veilleux & Rupke 
2002 1 . The line ratios harden with increasing height 



above the disk place. This could indicate the gas be- 
comes increasingly shocked at greater distance above the 
disk (not unexpected since th e wind spee d will increase 
with reduced ISM gas density (|Rand|1998[ )), or may sim- 
ply be a product of dilution of the outflow emission line 
components blended with the underlying disk emission in 
the simple single Gaussian line apporoximation currently 
implemented. The line ratio maps indicate harder line 
ratios on the western side of the southern lobe, coinci- 
dent with the x-ray emitting gas presented by , Strickland] 



NGC 253 

Due to its proximity, NGC 253 is one of the best stud- 
ied examples of the class of luminous infrared galaxies 
know n to be undergoing an intense burst of s tar forma- 
tion ([Rieke, Lebofs ky, Thompson aZ.|1980|). |Demouliii| 
|& Burbidge (1970) interpreted the disturbed kinemat- 
ics of the e mission line gas as a conic outflow which w as 
later shown ( McCarthy, van Breugel & Heckma n|1987[) to 
be_gualitatively consistent within the [Chevalier fe Ciegg 
(|1985| starburst wind model. Evidence for entrainment 
(rf mo lecular gas with in the outflow has been reported 
(Sakamoto etal. 2006). The outflow is clearly visible in 
the UV although the observed flux is most likely light 
from the cen tral starburst scatte red by dust entrained in 
the outflow ( Hoopes aL||20(j5 ). X-ray observations do 
not reveal evidence tor an AG N , being consistent with a 
starburst driven galactic wind ( Persic e^a/.|1998 Strick- 



land etal. 2000^ 

I'he large spatial extent of NGC 253 meant that only 
the southern side of the galactic plane could be surveyed 
with the SPIRAL IFU in the allocated observing time 
(^4 hours per object). The observations are presented in 
Figs.[5]-j7] NGC 253 closely follows the emission pattern 
of NGC 1482. Two prominent flngcrs of emission, in- 
terpreted as the limb brightened edges of a line-emitting 
cone, are clearly identified in the SPIRAL data (Figs.jsj), 
with the structure most clearly highlighted by the con- 
tours of [Nil] emission. Like M82, the [Oiii] emission is 
more centrally condensed than the low ionisation emis- 
sion (cf. Shopbell & Bland-Hawthorn 1998). Much like in 
NGC 1482, the prominence of the [Nil] emission is sug- 
gestive of shock excitation at the boundary of the flow. 

Line splitting, indicative of a hollow conic outfiow 
structure, is seen in spectra taken from the central region 
of the outfiow above the plane of the disk. The veloc- 
ity separation is roughly ~100 km s^^ close to the disk 
but rises to >230 km s~^ at the edge of the IFU mosaic 
suggesting an acceleration of the outflow with increasing 
height above the disk plane (under the assumption of a 
constant cone opening angle). The limb brightened cone 
edges show only a single velocity component. The wind 
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filaments converge to a considerably more compact re- 
gion ('^ 140 pc) in the nuclear disk compared to what is 
seen in NGC 1482. 



([2000[ ), the western 
In contrast, we find 



Consistent with Strickland etal 
limb is brighter in the Balmer fines, 
that the eastern limb is consistently brighter in all forbid- 
den lines relative to Ha (Fig.|6|. A similar asymmetry is 
seen in NGC 1482. We observe an apparent hardening of 
the line ratios with height above the disk, but this may 
be due to dilution from the inner disk. 

The IDDs constructed from the SPIRAL data (Fig. [7| 
reveal the dominance of shock ionisation in the optical 
filaments. In particular, weak [O i] emission is observed 
across the cone being especially prominent in the east- 
ern limb. At large distances, the weaker [Oi] emission 
is compromised by the proximity of strong telluric [O i] 
emission due to the low rcdshift of NGC 253. 

Consistent with this picture, Strickland etal. { 2000[ ) 
conclude that the outflow emission seen at x-ray energies 
is derived from a low filling factor gas that is excited by 
interaction between the wind and the ISM on the surface 
of the outflow cone, rather than from the wind itself. In 
conclusion, NGC 253 exhibits a starburst driven, shock 
excited wind nebula, observed in this instance only on 
the near side of the disk due to the spatial coverage of 
the data. The wind originates from a centrally conflned 
starburst, rather than from a signiflcant portion of the 
disk region. 

NGC 1808 

The "peculiar" nuc l ear sp ectrum of this object led 
Veron-Cetty fc Veron| ( |1985[ i to classify NGC 1 808 as a 
weak Seytert 2 galaxy Subsequent observations (Forbes, 
Boisson fc Ward|1992 ) showed that the properties ot the 
nuclear emission are more plausibly explained by a nu- 
clear starburst, a vie \ y supported by mid-IR diagnostics 



( [Laurent etaL|2000[ ). |Phillips| ( |1993[ ) confirm the nuclear 
outflow suspected from earlier work, arguing against the 
Seyfert 2 identification in favour of a starburst-driven 
outflow source. 



Junkes etal. (1995) flnd that the nuclear x-ray emis- 
sion from the unresolved core is largely due to the SNRs 
and hot bubbles of the central starburst. Near-IR spec- 



trosco py of the nucleus (Krabbe, Sternberg fc Genzel 
1994 1 favours a relatively young starburst (age ~ 50 iVlyr 
with a star formation rate of ^1 yr~^); see §4.1). 
But at higher energies, there is evidence for a heavily 



Jimenez-Bailon et al 
ert 2 classification in 



obscured AGN component (q.v, 
p)05)). Thus, we retain the Sey 
Table 2 for this galaxy, although the designation remains 

uncertain. 

Figs. [8l-[T0| present the new SPIRAL observations. 
NGC 1808 was observed under less than ideal conditions, 
with variable seeing and intermittent cloud cover. But 
we were able to extract useful IDDs for the inner extent 
of the wind. O ur observations reproduce the asymmetric 



line profiles of Phillips (1993), with blue and red shoul- 



ders corresponding to the near side outflow and counter 
flow, as previously reported. Due to the weakness of 
[O III] and [O i] , the IDDs do not provide an unambigu- 
ous answer for the nature of the ionizing source in the 
extended filaments. But we stress that the base of the 
wind is enhanced in [N ii]/Ha, similar to what is observed 



in NGC 1482, and the line profiles are somewhat broad- 
ened. These observations are consistent with large-scale 
shocks in the outflowing gas. 

NGC 3628 



Fabbiano, Heckman fc Keel 
strong evidence from x-ray anc 



(|1990[ ) presented the flrst 
Ha observations for a col- 
limated outflow suspected from earli er kinematic work 
(Schmelz, Baan & Haschick 1987a|b ). Subsequently, a 
cor respondmg mo l ecular disk- halo outflow was reported 
by llr^inj^ Sofue| d 1996[ ) . IR ( [Rice e^a/ |[l988[ ) and ra- 
dio ([Condon et a/.|1982p measurements place this object 
on tne well known tar-lR vs. radio correlation indicating 
that the central energetics are dominated by star forma- 
tion. At x-ray energies, no evidence is fou nd for an ac- 
tive n ucleus at the current detection limits ( Flohic et al 
2006| 



Ubservations of NGC 3628 were undertaken in less 
than ideal observing conditions. We were able to ex- 
tract useful spectral diagnostic information from the red 
data, but only limited information from the blue data. 
Emission line maps were extracted for Ha, [N ii] and [S ii] 
and are presented in Figs. [TT] & [12] The strongest line 
emission originates in a series of compact regions to the 
south-east of the SPIRAL mosaic, which we interpret as 
Hii regions. The [Nii]/Ha ratio image (Fig. 12) reveals 
the fllame ntary emission above the plane of tlie disk re- 
ported by Fabbiano, Heckman fc Keel (1990). 

The line ratios in this region show enhanced [Nil] /Ha 
similar to what is observed in most of our starburst wind 
sources. Strong [Oiii]/H/3 is ruled out by the blue data 
even with its compromised sensitivity. The base of the 
wind is enhanced in [Nil] /Ha, reminiscent of what is 
observed in NGC 1482, and the line profiles are some- 
what broader than line emission from neighbouring HII 
regions. Once again, these observations are consistent 
with large-scale shocks in the outflowing gas. 

NGC 6810 

The historic classification of this object as a Seyfert 2 
galaxy is based primarily on the width of emission lines 
in low spatial resolu tion spectra of the nuclear region. 
Coccato et al. ( 2004 ) noted "a higher than expected mi- 
nor axis velocity dispe rsion and unusua l kinematics" in 



long slit spectroscopy. Strickland (20071 concluded that 
NGC 6810 is in fact a "lull disk" superwind galaxy 
and that the previous classification is based largely on 
misidentification of counterparts in a low resolution ra- 
dio survey and the merged line profiles of an outflowing 
wind. This is co mpatible with the earlier observation 
( [Forbes fc Norris[ [T998) that the radio, FIR and [Fell] 
properties ot the inner disk are consistent with star for- 
mation. Our new SPIRAL observations (Figs. 
are consistent with this interpretation. There is no evi- 
dence for AGN activity in the optical spectra, with nar- 
row lines consistent with solar-metallicity star formation 
across the stellar disk. Star formation across the disk of 
the galaxy dominates the IDDs. 



The [O III] /Ha line ratio map (Fig. 14 1 shows strik- 
ing evidence for an "ionisation cone" on the near side 
of the disk. These structures are not unique to AGN 
sources (e.g. Pogge 1987; Tadhunter fc Tsvetanov 1989) 
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having been observed in starburst galaxies (e.g. Shop- 
bell & Bland-Hawthorn 1998). Leaked radiation from a 
starburst or an AGN, or even shocks excited by a conic 
outflow, can produce ionisation cones along the minor 
axis. In support of this statement, the [O ill] / H/3 ratio 



in th is reg ion is characteristically low (cf. jShuU fc Mc 
Kee (19791). In the line ratio map, we have used Ha in 



place of tip because the stronger Ha signal reveals the 
cone more clearly. 

Unfortunately, the weakness of [O ill] emission above 
the galactic disk makes it difficult to unambiguously clas- 
sify the excita tion mechanism for the outflow gas from 
the IDD (Fig. fis]). But we stress that the base of the 
wind is enhanced m [N ii]/Ha similar to what is observed 
in NGC 1482. The high [Nil] /Ha gas dehneates the 
boundary of a chimney structure emanating from the 
central starburst region. The base of the chimney ex- 
tends ~1.5 kpc across the inner galaxy. An inspection 
of emission-line profiles in regions of enhanced [Nil] /Ha 
reveals multi-component line profiles indicative of high 
velocity material in a blueshiftcd outflow on the east- 
ern side of the disk, with the corresponding counter flow 
redshifted (and dust attenuated) on the opposite side. 
In summary, our observations provide further evidence 
that the off-nuclear emission in NGC 6810 is due to a 
large-scale galactic wind. Like NGC 1482, the dominant 
source of ionisation is most likely due to shock excitation 
of filaments entrained in the wind. 

3.2. AGN 

NGC 1365 



Optical spectroscopy of this well known Seyfert dVeron, 
Lindblad, Zuiderwijk et at] |1980| [Edmunds & Page" 



1982 1 reveals a broad component to the Ha line and 
line ratios indicative of AGN activity. Further evidence 
for the AGN com es from x-ray varia bility of the highly 
obscured nucleus (Ward etal. 1978). But based on an 
analysis of t he radio and x-ray structure of the cen- 
tral regions, Stevens, Forbes & Norris| (|1999|) suggest 
that the energetics are dominated by star formation, and 
that "the evidence of a j et emanating from the nucleus 



is at best m arginal" (cf. Sandqvist, Joersaeter & Lind 
blad (1995)). The radio/ tar-lK ratio i s also indicative 



ot a d ommant circumnuclear starburst (|Forbes & Norris 
T998| 



The extended ionised gas in NGC 1365 has been exten- 
sively studied with long-slit spectroscopy . On the basis 
of these observations, ,Hjelm fc Lindblad ( 1996| ) propose 
a physical structure (wide angle, hollow conic outflow) 
to explain the unusua.1 v elocity structures previously re- 
ported (Phillips eiaZ.|1983; Joersaeter etal. 1984) in the 
context of an outflow (flrst suggested by Burbidgc & Bur-' 
bidge (I960)), while _Storchi-Bergmann fc Bonatto^(1991i ) 
present clear evidence ot an [O iiij "ionisation cone" to 
the south east of the nucleus, above the plane of the disk. 

The new SPIRAL observations for NGC 1365 are pre- 
sented in Figs.|16l- )18[ The reconstructed continuum im- 
age reproduces the central disk structure seen in the HST 
imaging from Fig. [l] The Ha integrated flux map sug- 
gests strong star formation in the nuclear zone, with a 
number of probable compact Hll regions or star clusters. 

Strong [O III] emission is seen to extend eastwards from 
the nuclear region above the plane of the disk, while a 



counterpart emission region seen below the disk is at 
lower intensity, presumably due to dust attenuation on 
crossing the plane of the nuclear disk. [N li] emission ex- 
hibits the classic limb-brightened arms of a conic struc- 
ture to the west but such a feature is less obvious on 
the near side of the disk, possibly d ue to confusion with 
disk emission. Phillips etal. (1983) report the "remark- 
able kinematics ot the ionised gas m the nucleus of NGC 
1365" and line splitting to be most prominent in the 
south-eastern region, particularly in the [O ill] line, ob- 
servations supported by the SPIRAL data. Line splitting 
is evident but not complete at the resolution of the SPI- 
RAL observations, suggesting a narrow opening angle or 
low velocities for the outflow. 



The emission-line ratio maps shown in Fig. 17 
revealing. The counter fan of [O ill] emission is evi- 
dent to the west in the [Oiii]/Ha ratio, while the en- 
hanced [Nil] /Ha ratio usually observed in outflowing 
emission line gas is seen to flll the region bounded by 
the limb brightened [N ii] structure of the western region 
of Fig. 16 The expected counter-flow region is visible to 
the east as high [Nil] /Ha gas once far enough above the 
disk to avoid contamination by the strong star formation 
in the region. The [S ii] electron density diagnostic shows 
no clear structure at the limit of our signal-to- noise limit. 
The ionisation cone from the central AGN is most evi- 
dent on the nearside of the disk in the [O ill] / [O i] ratio 
{Fig.m. 

The classic IDDs are presented in Fig. [18} the dom- 
inant signature is that of a high metallicity starburst 
consistent with dominant circumnuclear star formation. 
However, significant emission, particularly in the outflow 
re gions, lies above the extreme star formation boundary 
of 
al 

of the hard ionising spectrum expected from AGN pho- 
toionisation, and is for the most part inconsistent with 
photoionisation by starlight or shock ionisation models. 
Broad Heii A4686 and Hei A5876 emission are observed 
in the nuclear spectrum, consistent with photoionisation 
by a hard ionizing source. 



CIRCINUS 

Located at low galactic latitude and with high fore- 
ground extinction (Ay =1.5), t his large nearby Seyfert 
galaxy was first identified by Freeman et al.\ |i977). 
Maser measurements firmly est ablish the central" black 
hole mass at ttibh ~ 2 x 10^ Mq ( jCreenhifl al.\2003al 
With the inferred accretion rate (ftiU.l Lpdd ) at odds with 



Kewley etal. (2001 ). The zone occupied by this gas, in 
three of the classic diagnostic diagrams, is indicative 



:!U.l Ledd) at odds witl 
much higher accretion rates in other Seyferts with maser- 
determined black hole masses (e.g. NGC 1068; Begelman 
fc Bland-Hawthorn 1997), Circinus lends support to the 
notion that supermassive black holes are indeed quite 
common in galaxies, but for the most part have long 
duty cycles and/or low accretion rates. 
The Seyfert 2 nature of Circi nus is supported b y the 



mid-IR coronal li ne emission ([Lutz et al. 2002), e.g. 
[Siix] at 3.94/zm ([Oliva etaI1|l994p. Recent integral 



field spectroscopy, assisted by adaptive optics (Mueller 



Sanchez, Davies, Eisenhauer et al. 2006 ), indicates a sub- 
stantial (1— 2%Lboi) young (<1(J0 Myr) toroidal star- 
burst within 10 pc of the central AGN and provides 
further evidence for a "wide angle outflow" as reported 
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by Greenhill etal. (2003b I. Extended filamentation as- 



sociated with the~buttlow was first mapped by IVeilleux 
fc Bland-H awthorn] (11997 ) , and later in more detail by 
Veilleux et al. (2003^ X-ray spectral analysis further 
supports the view oF^Circinus as a highly obscured AGN 
but with strong thermal dust emis sion dom inated by 
heating from nuclear star formation (Matt etal. 2000). 

The large spatial extent of the Circmus galaxy allowed 
only the north-west side of the disk to be probed with 
SPIRAL. Figs . [l9l-[2T] present the observations. Our spa- 
tial covera ge was guiaed by TTF observations of |VeilTeux 
fc Bland-H awthorn ( 1997) and the HST study of |Wilson, 
Shopbell~Simpson^ (2000j). We have not atteinpted to 
correct the observed spectra for for eground reddening 
The "hook" stru cture identified by " ~' 



Veilleux & Bland- 



Hawthorn ( 1997 ) traces the southern boundary of an ion 
isation coiie, running east-west in the prominent emission 
lines. 

The [O III] line profiles show significant complexity 
with a number of mixed velocity components. Maximal 
line splitting in the centre of the ionisation cone, and sin- 
gle line profiles along its outer walls, might suggest a sim- 
ple conic outflow structure. But the line intensity maps 
indicate limb darkening rather than the limb brightening 
one would expect on a conic emitting surface at the inter- 
face of a hot wind material and the cooler ISM, suggest- 
ing significant emission from gas within the cone. There 
is a strong velocity gradient across the filaments, each fil- 
ament being at a different mean velocity, with the [O ill] 
emission from the east-west boundary significantly red- 
shifted with respect to the opposite (north-south) region 
of the ionisation cone. The onflow is best interpreted 
as a gas-filled ionisation cone illuminated by the central 
AGN, rather than a walled cavity structure. 

NGC 5128 (CENTAUR US A) 

Cen A is one of the most well studied active galaxies 
owing to its extraordinary brightness across 16 decades 
or more of the electromagnetic spectrum (Israel 1990). It 
is the nearest powerful AGN and, more specifically, the 
closest FR I radio source. Fortuitously, we are witnessing 
an accretion-powered radio source in the first 100 Myr of 
its projected 0.5 — 1 Gyr lifetime. A well known feature 
of Cen A, unique amo ng our sample, is the presence of a 
spectacul ar radio jet ( Blanco, Graham, Lasker &: Osmer 
[l975 , Peterson, Dickens & Canno n II197 5[) which ignites" 
optical filame nts along its axis ('Sutherland, Bicknell & 
Dopita |1993 ). The jet- induced emission falls tar beyond 



the radial extent of our IFU mosaic. 

Over the inner galax y, the ex tended ionised gas was 
first mapped by Bland, Taylor fc Atherton (1987) and 



Nicholson, Bland-Hawthorn & I'aylor (1992|, and shown 
to arise exclusively trom a severel y-warped thin disk 



first : 



later confirm ed by CO observations (|Quillen, Graham & 
Frogel||1993|) 
m 



A diffuse highly ionised halo of [O liij was 
entified by Bland-Hawthorn & Kedziora-Chudczer 
(2003, see their Fig. 6) and has recently been sho wn to 
match the extended x-ray maps from Chandra ( Kraft 
etal. 2008) extraordinarily closely. Extended niid-lR 



emission is also seen along the rad io axis, revealing ev - 
idence for a hmb-brightened shell (Quillen etal. 2006), 
and providing strong evidence for a bipolar outflow along 
the radio axis. 



Here we provide only a brief summar y of the SPIRAL 
observations presented in Figs. [22[-[24j A more detailed 
discussion and associated modeling is to be presented 
elsewhere. The [Oiii]/H/3 vs. [Nii]/Ha IDD reveals a 
large number of data points above the extreme starburst 
line. When these points are identified in the SPIRAL 
data, they reveal a striking nuclear ionisation cone. The 
cones vissibility is enhanced after a 3x 3 spatial binning of 
the mosaiced blue IFU data (Fig. 25 ) . The conic region, 
delineated by sharp edges, exhibits gas with the ioni- 
sation ratios expected from a hard AGN ionising spec- 
trum, although shocks cannot be ruled out. The struc- 
ture can be traced back towards the nucleus wh ere it 
likely connects with the Pa/? em ission cone seen by |Kra- 
jnovic. Sharp & Thatte (2007) at near-IR wavelengths. 
I'he feature is roughly aligned with the known radio jet 
structure above the plane of the disk. Our coverage to 
the south is insufficient to trace a possible counter cone. 
Throughout the cone, we clearly see variations in the 
[S ii] doublet consistent with density variations of order 
rig PS 300 — 1000 cm~'^. This is consistent with the in- 
ferred shell -enhanced den sities from the mid-IR observa- 
tions ( Quillen a/.||2006[) . 

In summary, we have discovered an "ionisation cone" 
along the minor axis of Centaurus A for the first time. 
The phenomenon is most evident in emission-line ratio 
images (e.g. [Oiii]/H/3) where a clear distinction can be 
made between the ionising source in the outfiow from 
the background ionising source. This tends to favour 
an ionisation cone powered by an AGN source, but the 
phenomenon can also be identified in starburst galaxies 
(see §4.1), as we discovered here in NGC 6810. 

IC 5063 



Bergeron, Durret & Boksenberg ( 1983 ) presented early 
evidence tor an obscured Seyfert nucleus via the detec- 
tion of a broad yet weak component to the Ha emis- 
sion. iColina, Sparks & Macchetto ( 1991 ) reported high- 



excitation emission lines and confirmed the broad Ha 
profile, thereby supporting the claim of a concealed ac- 
tive nucleus. Both papers note the kinematic decoupling 
of gas in the inner regions from the extended disk rota- 
tion c urve of the galaxy. |Colina, Sparks fc Macchetto 
( 1991 1 favoured an explanation invoking outflowing gas, 
and provided evidence of a wide angle (50°) ionisation 
cone via the spatial map of the [O ill] /Ha ratio . 

At radio wavelengths, Morganti etal. (20071 present 
high quality observations oFwarm and c ool material mov- 
ing outwards from the nu cleus (see also Morganti, Oost- 
erloo fc Tsvetanov (1998)). Their observations suggest a 
radio lobe/jet structure which in turn may be responsible 
for the disturbed kinematics of the outflow gas. These 
authors conclude there is little evidence for shock-driven 
excitation, proposing instead that the emission is over- 
whelmed by AGN photoionisation. 

The SPIRAL observations a re shown in Figs. [26 
As presented in previous work ( Bergeron, Durret fc 



senberg 1983 Colina, Sparks "fc" Macchetto 1991j), the 



velocity field ot the emission line gas is dominated by 
a regular disk rotation beyond 2" from the nucleus, al- 
though line profiles show evidence of broadening at the 
base of the stronger lines, and line splitting is apparent 
above the mid-plane. The kinematically distinct outfiow- 
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ing gas shows little evidence of the classic line splitting 
expected from a hollow conic outflow structure, indica- 
tive of a more homogeneous gas distribution in this re- 
gion. The ionisation diagnostic diagrams (Fig. |28[ ) in- 
dicate AGN photoionisation almost exclusively, at least 
in regions where the intensity of emission is sufficently 
above the fitting threshold. 

NGC 6240 

As one of the closest of the ultraluminous infrared 
galaxies (ULIRG), NGC 6240 is central to the longstand- 
ing debate as to the presence or absence of quasar-like 
nuclei in ULIRG sources. In the case of NGC 6240, hard 
x-ray observations unambiguously identif y a double AGN 



core at its heart (Komossa etal. 20031. A large-scale 



wind was first reported by Heckman, Armus and Mi 



ley (1990) (cf. Ptak, Heckman, Levenson etal. ( 2003 ), 
Colbert, Wilson Bland-H awthorn' (1994|). The ex- 
tended emission-fine gas has been mappe d m optical lines 
by iniaging tunable filters ([Bland-Hawthorn, Wilson fc^ 
TuIIylfiggil IVeilleux etaL||2003[ ) and in x-ray emission 
( [Koniossa eit a/.||2003p . 

I'he SPIRAL observations show that the nuclear re- 
gions of NGC 6240 are highly complex both in the kine- 
matics and ionised gas distribution. With limited spatial 
coverage possible due to observing c onstraints, our obser- 
vations were guided by the study of Veilleux et al. ( 2003 1 
and focus on the region of high [N Iij/Ha ratio eniission. 
Two dimensional emission maps are not shown as the 
complexity of the emission structure in the nuclear region 
combined with the low surface brightness of emission in 
the less complex outer zones is not easily interpreted in 
this manner. Spectra from the nucleus show kincmati- 
cally broadened line emission (FWHM-13A) with a ve- 
locity spread of 500-600 km s^^. Most spectra do not 
easily separate into identifiable Ha and [N ii] lines. 

An emission-line filament is observed to extend east 
from the nucleus. We interpret the isolated line in the 
filament spectrum as Ha offset in velocity from the host 
galaxy by -1-340 km s~^in the outer reaches of the mosaic. 
This is at odds with the Ha vs. [N ii] map of Veilleux et al. 
(2003) but we note that the proposed shift in the emis- 
sion wavelength could lead to cross-contamination be- 
tween the wavelength settings of the tunable filter obser- 
vations. The conspicuous absence of satellite Ha (and 
the fainter [N ii] line) in the filament spectrum support 
our proposed identification. On heavily binning adjacent 
spectra, there is evidence for a weak feature compati- 
ble with broad [Nil] emission to accompany Ha. There 
is tentative evidence for [O ill] emission associated with 
the Ha filament as well, although the sensitivity of the 
data is not sufficient to confirm this in detail. In sum- 
mary, we are unable to construct IDDs for this object 
with the present observations. 

4. AGN VS. STARBURST 

In the previous section, we saw consistent differences 
in the ionisation characteristics of outfiowing gas associ- 
ated with starburst and AGN sources. We now present 
a simple model that largely explains these differences. 
Consequently, our model provides a key insight into the 
galactic wind mechanism and the radiative timescales of 
the starburst and AGN phenomenon. 



4.1. Predicted and observed levels of photoionisation 

Our galaxy sample spans bolometric luminosities in the 
range Lboi = 1 — 10 x 10^° L© and are distinguished by 
their nuclear starburst or AGN activity. All of the galax- 
ies were selected on the basis of known minor-axis out- 
flows. The new observations reveal that all of the AGN 
galaxies display IDDs typical of direct photoionisation by 
the AGN, whereas the starburst galaxies display IDDs 
typical of shock ionisation. How are we to understand 
this? 

First we note that both nuclear starburst and active 
galaxies almost certainly have duty cycles such that what 
we observe now is more representative of their peak activ- 
ity rather than their quiescent phase. We now show that 
prima facie the distinction between starburst and AGN 
outflow IDDs is surprising when one considers the overall 
ionizing photon production in both classes of sources. 

A black hole that converts rest-mass energy with an 
efficiency e into radiation has a luminosity 



La — erne 

V0.05/ 



Mq yr-1 



Lr. 



(1) 

(2) 



for which m is the mass accretion rate. In reality, the 
accretion disk luminosity can limit the accretion rate 
through radiation pressure. The so-called Eddington 
limit is given by 



Le = 



AnGmBii'mpC 



= 2 X 10 
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™bh 

W Mr: 



Lr. 



(3) 
(4) 



where mBH is the black-hole mass, rup is the proton mass 
and ax is the Thomson cross-section for electron scat- 
tering. This equation allows us to derive a maximum 
accretion rate to, and therefore a maximum ionizing lu- 
minosity, for a given black-hole mass. We find that this 
relation is not strongly constrained with the present data. 
For a canonical power-law continuum (oc {hu)~^) typical 
of active galactic nuclei, the number of Lyman continuum 
(LyC) photons is given by 



LyC,A 



Vio"i, 



phot s 



(5) 



for which is the fraction of Lboi {La — S,aL\,o\) aris- 
ing from AGN activity today (i.e. not averaged over the 
lifetime of the source). 

The total radiative output from a luminous starburst 
at the present epoch can be approximated by (e.g. Mal- 
oney 1999; Heckman, Lehnert & Armus 1993) 
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5 Mq yr-i 

0.23 



/ At 



1 M, 







771(7 



(lOO Mq 



\ 108 yr 

0.37 



0.67 



Lr. 



(6) 
(7) 



where Ai is the starburst duration estimated from the 
gas supply divided by the star formation rate s. The 
third and fourth terms indicate the lower and upper stel- 
lar mass cut-offs. For a conventional IMF, the star for- 
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mation rate is roughly 



(8) 



We can now determine the rate of Lyman continuum 
photons produced by the starburst. To a useful approx- 
imation, 



LyC.S ' 



10 



54 



5 Mq yr- 

Lhoi 
IQIIL; 



rriL 
1 AU 



phot s ^ 



(,100 Mo j 

(10) 



which is to be compared to equation [5] above. The high- 
energy cut-off for the hot young stars declines rapidly 
beyond 50 eV. The factor ^5 is the fraction of Lboi {Ls = 
Cs^boi) arising from the starburst activity today (i.e. not 
averaged over the lifetime of the source). 

The LyC photon rate for both starbursts and AGNs 
for two sources with the same bolometric luminosity is 
essentially the same. While the nonthermal spectrum of 
the AGN means that the harder photons propagate at 
least an order of magnitude further into the surrounding 
gas, the photon number is dominated by the soft photons 
so that the comparison is still meaningful. We do not dis- 
cuss the possibility of the total accretion disk luminosity 
being channelled into a smaller solid angle (i.e. intrinsi- 
cally beamed) because we were unable to find compelling 
evidence for this to be so. This is to be addressed in our 
later papers. 

An ionizing photon rate as high as AfhyC ^ 10^"*^ 
(Lboi/lO^^i©) phot s^^ is easily large enough to account 
for the off- nuclear line fluxes observed in all galaxies, 
whether starburst or AGN. At a distance of 1 kpc, the 
ionizing flux is 1^9 10^" phot cm~^ s~^. We can relate 
this to the Ha surface brightness through the emission 
measure where 



f 



fn^ dz cm pc 



(11) 



which is an integral of H recombinations along the line 
of sight z multiplied by a filling factor /j^ If we assume 
that the absorbing cloud along the minor axis is opti- 
cally thick, we can relate this directly to an Ha surface 
brightness in cgs units such that 



^(Ha) = 1 X IQ- 



1 kpc 



L 



bol 



lOiiL. 



Cgs (12) 



where r is the distance to the central source (Bland- 
Hawthorn 1998, Appendix). Therefore the expected 
Ha surface brightness at 1 kpc is roughly 7 x 10~^^ cgs 
for the AGNs, and 2 x 10~^^ cgs for the starbursts. (The 
difference here relates to the different mean bolometric 
luminosities between the two classes of objects.) These 
predictions are an order of magnitude brighter than what 
is seen in our objects (Table 4), with the possible excep- 
tion of the AGN IC 5063 where the minor axis filaments 
are very bright. The higher predicted values may be due 

^ £m(Ha) = 1 cm~® pc is equivalent to 2xl0~^** erg cm~^ 
arcsec"'^ (cgs) for a plasma at 10'' K. 



to a combination of several factors: (i) the emitting re- 
gions have a small filling factor (/ < 1), (ii) some dust 
is mixed in with the filaments (cf. Cecil etal. 2001), and 
(iii) ^ < 1 for the AGN/starburst fraction of the bolo- 
metric luminosity. 

Thus we find that both AGNs and starbursts have the 
capability to ionise much of the extended gas filaments 
observed far along the minor axis, assuming that the ra- 
diation is not scattered or absorbed en route. So we ar- 
rive at an important conclusion. The distinction between 
starburst and AGN outflow IDDs is one of timescale, as- 
suming that the wind mechanism is not too dissimilar 
between the two sources (e.g. VeiWeux. et al. 2005). The 
time evolution of UV photon production must be radi- 
cally different in both cases. To explore this issue further, 
we examine more closely how the energy is deposited into 
the surrounding medium. 

4.2. Time evolution of UV and mechanical luminosity 

in a starburst 

The starburst phase in a galaxy's life is a spectacular 
phenomenon that can be observed to very high redshifts. 
Just how the formation of hundreds of millions of stars 
can be triggered over the central regions of galaxies in 
a dynamical time remains an open issue. What is well 
known is that the high mass stars evolve rapidly, giving 
off prodigious amounts of UV radiation initially before 
entering into the supernova phase, whereupon a huge 
amount of mechanical energy is deposited into the ISM 
(Leitherer etal. 1992). 

In Fig. |29[ we show the results from running Star- 
burst99 with a conventional initial mass function (IMF) 
for both instantaneous and continuous star formation 
models. We show the time evolution of the ionizing and 
mechanical luminosities arising from the starburst. Two 
upper mass cut-offs are presented in order to illustrate 
the imp ortant contribution of the highest mass stars. In 
Fig. 29 ^a) , the rapid drop in UV production at 10 Myr 
is evident, as is the rise of mechanical output from su- 
pernovae beyond here. This corresponds to stars more 
massive than IZMq evolving off the main sequence to 
b eco me s up ernovae . 

Fig. 29 a) can readily account for why stellar photoion- 
isation does not dominate over shock ionisation. Thus, 
we conclude that the ionisation timescale (T/) ind icat ed 
by starburst winds is such that T/ > 10 Myr. Fig. [29| b) 
effectively rules out the continuous star formation rnodel, 
at least for the starbursts in our sample. The model pre- 
dicts that stellar photoionisation is detectable at all times 
contrary to what we observe. 

In the next section, we show that the predicted me- 
chanical luminosity in the instantaneous star formation 
model is more than enough to explain t he o bserved shock 
ionisation. Indeed, we argue that Fig. [29] dates the ob- 
served starburst wind phenomenon in the range 10 to 40 
Myr. A more precise age will depend on how the star- 
burst region evolved early in its life (Suchkov et al. 1994; 
Cooper etal. 2008). We discuss this issue in more detail 
below. 

4.3. Time evolution of UV and mechanical luminosity 

in an AGN 

We must now explain the dominance of AGN photoion- 
isation of the filaments in AGN wind objects. There 
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are few reliable constraints on the long-term temporal 
evolution of accretion-disk radiation fields. Like star- 
bursts, AGNs clearly display all the hallmarks of duty 
cycles, bursts and outflows over the course of their life- 
times (Soltan 1983; Cavaliere & Padovani 1990). Most 
galaxies appear to have central massive black holes, but 
only a few percent can be classified as AGNs. The most 
likely explanation is that almost any galaxy can display 
AGN activity if its nuclear regions experience substantial 
gas accretion from outside. 

But we arrive at the long-standing debate on what 
can be uniquely attributed to the AGN as compared to, 
say, circumnuclear star-forming activity (Terlevich et al. 
1985). After all, the gas inflow that triggers the former 
can just as easily give rise to the latter. These activities 
are commonly seen in close proximity, in particular, in 
Seyfert galaxies. We address both oiF these issues now, 
i.e. why the AGN dominates the photoionisation of the 
wind filaments, and whether our conclusions are affected 
by the presence of the circumnuclear starburst. 

A reasonable AGN scenario is that the activity is ini- 
tially triggered by a gas dump, followed by a slow decline 
in the fuelling process (e.g. Norman & Silk 1983) as the 
gas settles down within the nuclear potential well. We 
assume an exponential model for the UV output from an 
AGN such that 



La/L^oX = C'A^maxexp( — t/t^) 



(13) 



where CA,max is the maximum fraction of the bolometric 
luminosity that the AGN achieves over the lifetime of the 
activity with half life t^. We further restrict this relation 
using the following normalization: 



LL 



La 
hv 



dv = TVlyCA 



(14) 



where LL denotes the Lyman limit (hv ~ 13.6 eV) and 
the RHS is given by equation [Sj In practice, this normal- 
ization for CA,max and tA is oifRcult to do properly be- 
cause our observations constitute an instantaneous snap- 
shot, rather than the time- averaged behavior of its prop- 
erties. Of course, the same holds true for starburst activ- 
ity. In principle, IFU studies of a large sample of AGNs 
and starbursts could provide accurate measurements of 
£,A and and therefore provide strong constraints on 
C'max for both AGNs and starbursts. 
If we integrate equation 13 over time (At ~ 100 Myr), 



subject to the normalization m equation |14[ we arrive at 



the important result that Ca,] 



cx t. 



The normal- 



ization ensures that a range of AGN models depending 
on CA,max and Ia produce the same ionizing radiation 



field averaged over the lifetime of the source. In Fig. 30 



we superimpose four s uch models on the instantaneous 
starburst model in Fig. [29j 

As already mentioned7 essentially all of our AGN 
sources show appreciable circumnuclear star formation. 
Since AGN photoionisation is clearly dominant over stel- 
lar photoionisation in all of our AGN sources, the pre- 
dicted surface brightness values effectively rule out AGN 
exponential timescales below 10 Myr. If the AGN UV 
radiation field is "beamed" along the wind axis, in the 
sense that the AGN source is strongly anisotropic, then 
shorter AGN timescales are possible. But the dynamical 
timescale of the wind discussed in the next section argues 



against this interpretation. Furthermore, we have been 
unable to find any reliable evidence for a UV boost along 
the wind axis in AGNs. 

But a question remains. In Fig. |30[ note how the 
starburst mechanical luminosity can dominate over the 
AGN UV luminosity beyond 10 Myr, assuming the 
mass/energy input can be efficiently converted to ion- 
izing radiation. Why do we therefore not see the signa- 
tures of shocks in AGN-driven winds? Shocks are seen in 
starburst winds because the central ionizing source fades 
rapidly early in the evolution of the wind. But this prob- 
lem only arises if we assume that the mechanical energy 
in the AGN sources is supplied with a long delay as for 
starburst-driven winds because, for example, the AGN 
wind is assisted by a circumnuclear starburst. The ob- 
servations may therefore indicate that the AGN-driven 
winds are fundamentally different from starburst-driven 
winds. In this scenario, the mechanical energy is sup- 
plied by the accretion disk immediately, but the shock 
ionizing luminosity never exceeds the luminosity of the 
central accretion disk. 

We now consider the shock ionisation observed in 
starburst-driven winds. We find that only a fraction 
of the available mechanical luminosity must be used in 
shock-exciting the gas, either because the process is in- 
efficient or because the filling factor of the extended fila- 
ments is lower than unity. 

4.4. Dynamical timescale of a large-scale wind 

Here we address the dynamical timescale of galactic 
winds before going onto derive levels of optical emission 
powered by shock ionisation. Once the starburst or AGN 
fires up, it begins to deposit energy into its surroundings. 
This largely mechanical energy over-pressurizes a circum- 
nuclear cavity of hot gas that reaches a temperature 



T « 3 X 10* X A" 



K, 



(15) 



where x is the efficiency of turning mechanical energy 
into heat, and A is the mass- loading factor (Veilleux et al. 
2005). This ratio measures the total mass of heated gas 
to the mass that is directly ejected by SNe and stellar 
winds or by the AGN. Most galactic winds appear to 
entrain gas from the interstellar medium as sho wn by the 
rotat i on of this material about the flow axis (see Veilleux 
et al. ( 2005 ) for a detailed discussion) . The entrained gas 
winch can easily dominate the wind mass, is accounted 
for with A typically larger than unity (e.g. Suchkov et al. 
1994). 

The central cavity is too hot to cool efficiently and so 
builds up a pressure 



Po/k - 2 X 10*^ 



5 Mq yr-iy Vl kpc 



K cm 



(16) 

where i?* is the radius of the energy injection zone. 
This pressure can significantly exceed the pressure of 
the undisturbed ISM. The hot cavity then evolves like 
a stellar wind-blown bubble (e.g.. Castor, McCray, & 
Weaver 1975; Koo & McKee 1992a,b). As hot gas ex- 
pands through the sonic radius, it cools adiabatically. 
Beyond the sonic radius, the wind drives a shock into 
the surrounding ISM and starts to sweep up a shell of 
shocked gas, a process that gradually slows the bubble 
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to much less than the wind velocity. This marks the 
end of the adiabatic "free expansion" phase, whose du- 
ration is set by the mechanical luminosity of the star- 
burst or AGN and the original density of the ISM. After 
free expansion, the system develops distinct concentric 
zones from the centre moving outwards due to the strong 
shocks (Cooper et al 2008; Dopita & Sutherland 1996; 
Weaver eiaZ. 1977). 

Energy-conserving bubbles. If radiative losses of the 
overall system are negligible, the expanding bubble is 
energy-conserving. In that case, the radius and the ve- 
locity of the expanding shell of shocked ISM are given by 
(Castor eiaL 1975; Weaver ei a/. 1977) 

rshcu = 1.1 (^44/«o)'/' 4^' kpc, (17) 

and therefore 

14hcii-640 (^44K)'/' i^'/' (18) 

= 670 {E,,/n,y/^ r;Zi'^^^ km s'^ (19) 

where is the age of the bubble in Myr, rio is the am- 
bient density in cm^"^. By convention, we use £'44 for 
the mechanical luminosity of the wind in units of 10'*'* 
erg s-i {Ei4 « x lO^Lo). 

Momentum-conserving bubbles. If radiative losses are 
significant, momentum conserving bubbles decelerate 
somewhat faster such that rghoii c>c t^;Vshcii (x (x 
''shell (s-S-) Avedisova 1972; Steigman, Strittmatter, & 
WiUiams 1975; Koo & McKee 1992a,b). For the galax- 
ies we are considering here, the inferred timescales are 
roughly comparable (see Table [5]). 

TABLE 5 

A COMPARISON OF EXPANSION TIME (tg IN Myr) AND SHELL 
VELOCITY (\4hell IN KM S"'^) FOR A BUBBLE THAT HAS EXPANDED 
TO A RADIUS rshell (iN KPC). VALUES ARE SHOWN FOR AN ENERGY 
CONSERVING AND A MOMENTUM CONSERVING BUBBLE. 

energy conserving momentum conserving 
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616 
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419 


4.4 


503 


4 
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346 
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435 
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8.1 


298 


12.2 


390 



In most simulations to date, the terminal velocity of 
gas entrained by a galactic wind is Vt ^ 400 — 800 km 
s~* (e.g. Suchkoveta^. 1994; Cooper a/. 2008). Proper 
comparisons are difficult because the models deal with 
dumpiness in the entrained gas, which the scaling rela- 
tions above do not. Chevalier & Clegg (1985) provide a 
useful formula for the terminal velocity of a cloud accel- 
erated by a wind (see also Strickland & Heckman 2008): 

V, « 430 kpc)-V2(aJ10-3g cm-^)-*/^ 

(20) 

where Tq and CTq are the initial distance and column den- 
sity of the cloud. (This formula is specific to the "con- 
stant velocity" phase of the wind before it breaks free of 
the confining medium.) The factors fe and fm account 



for the fractions of supernova energy and mass respec- 
tively that work to heat the diffuse gas. These authors 
explore fractions in the range 0.1 to 1 and ultimately 
favour the high end of this range. In any event, the dy- 
namical time of the wind is roughly tjj > ro/Vt ^ 2 Myr 
(see also Cooper et al. 2008) which is considerably longer 
than the cooling time of the entrained gas. 

The wind timescale is comparable to the dynamical 
time of the starburst region itself, but is up to an order 
of magnitude less than the ionisation timescale Tj estab- 
lished by our analysis above. This is a striking result that 
suggests a substantial delay {> 10 Myr) after the onset 
of star formation before the conditions for a large-scale 
wind are properly established. 

4.5. Shock ionisation in a large-scale wind 

In the starburst wind objects, the overall ionisation 
properties are consistent with shock processes (see §3.1). 
A simple calculation can validate this. From Dopita & 
Sutherland (1996), the total shock luminosity is 

Lshock « 2 X 10^*V;?,o^k'^o^hock erg s"*. (21) 

These authors find that about 50% of the mechanical 
energy is converted into energy capable of ionizing gas. 
Assuming there is sufficient gas to soak up the ionizing 
photons, one can derive a useful relation for the Ha lu- 
minosity, such that 

ishock(Ha) « O.OlLshock't^ehock erg s"*. (22) 

We now derive a canonical Ha surface brightness for 
shock-ionised filaments associated with the galactic wind. 
For a galaxy at a distance of £> = 10 Mpc, the observed 
Ha flux from the shock is 

/shock(Ha)«2.4 X 10-2V,L'ck"oAshock (^YoJa^) ^^^^ 

erg cm"^ s"* (24) 

which assumes, consistent with Dopita & Sutherland 
(1996), that only a single shock surface is observed with 
no correction for possible limb brightening. 

The surface brightness of shocked emission at 1 kpc 
follows: 

«„.,H„,. 5.2x10-.= (j5^)"(j^) cgs 

The pre-shock density normalization assumes that the 
halo gas in starburst galaxies follows an isothermal dis- 
tribution (e.g. Cooper et al. 2008) normalised to a pre- 
shock density of Uq ~ 10 cm^"^ at 1 kpc consistent with 
material swept up in a wind (q.v. Chevalier & Clegg 
1985; see below). Post-shock densities inferred from the 
[SII] doublet (e.g. McCarthy ei aZ. 1987) are typically in 
the low density limit (< 100 cm~'^) beyond r = 300 pc. 

While equation [25] has a strong dependence on Vghock: 
the shock speed is well constrained by observation. Some 
authors have associated the terminal velocity of en- 
trained gas with the shock velocity induced by the wind 
(e.g. Maloney 1999). But the levels of ionisation tell a 
different story. The strong [NII]/[OIII] signal favours 
shock velocities of order 100 km s~* (e.g. Shopbell & 
Bland-Hawthorn 1998; Veilleux & Rupke 2002). This is 
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supported by the soft x-ray emission which appears to 
be powered by similar shock velocities. How are we to 
understand this? 

In our picture, the fast hot wind accelerates denser gas 
to a terminal velocity. These clouds are ram-pressure 
confined by the fast wind. We therefore write down a 
simple equation for shocks driven into the clouds such 
that 



hock 



(26) 



where Uyj and are the wind density and velocity. 

Strickland & Heckman (2009) have recently conducted 
an excellent re-evaluation of Chevalier & Clegg's original 
model for the M82 wind. The free-flowing or break-out 
wind speeds can be as high as = 3000 km s~^. For 
our adopted values of Uo, Vshock and Vt above, we can 
balance equation 26 with ^ 0.01 cm~'^, a factor of a 
few higher than the canonical solution at 1 kpc presented 
by Strickland & Heckman (2009). We consider this to be 
broadly consistent given the uncertainties. In our simple 
analysis, the shock velocity is always several times less 
than the terminal velocity of clouds caught in the flow. 

Our simple model leads to an Ha surface brightness 
of w 5 X 10^"'^^ cgs at 1 kpc. This approximate value 
is comparable to what is presented in Table |4] for the 
starburst galaxies. In principle, the filling factor of the 
shocked material can be low (/ < 1), and there will be 
moderate levels of extinction internal to the entrained gas 
(i.e. Av < 0.3 mag; Cecil etaL 2001). But in principle, 
a moderately fast shock could account for the observed, 
extended line emission. 

We note that the Ha surface brightness in equation 25 
is comparable to the observed AGN surface brightness in 
Table |4j but much fainter than the predicted AGN sur- 
face brightness in equation |12[ This may indicate that 
the filling factor in the windfilaments is / ~ 0.1. For- 
tunately, a discussion of the relative merits of different 
ionisation sources does not require knowledge of the gas 
filling factor /. A lower filling factor affects the theo- 
retical prediction identically in equations [12] and |25[ and 
usefully brings our models into line with the observa- 
tions, at least in principle. Thus, if AGN winds involve 
similar processes to starburst winds, we would not ex- 
pect shocks to outshine direct photoionisation from the 
AGN, consistent with our observations. 

We have already argued that shocks outshine direct 
stellar photoionisation from the starburst because the 
young starburst has dimmed tremendously by the time 
the wind breaks out. Close to the disk, we may still 
expect stellar photoionisation to dominate over shocks 
owing to large-scale star formation in the disk (Bland- 
Hawthorn & Maloney 1999). But much beyond 300 pc 
or so, wind-driven shocks can dominate the ionisation, 
as observed in our sample. If we invoke a lower filling 
factor consistent with AGN winds, we can readily explain 
NGC 3628, but now NGC 253 and NGC 1482 appear 
overluminous. Indeed, a careful study of Table [4] reveals 
the need to change / over the range (0.1,1) to match our 
model to all of the observations. This may indicate that 
our model is oversimplified and needs to include a gas 
filling factor, or even a correction for internal extinction 
(although this is found to be typically less than 10% from 
the SPIRAL observations in the brightest regions). We 
note two things: (i) the inferred filling factors are not 



unreasonable and can be directly measured in principle; 
(ii) alternatively, the inferred extinction corrections are 
plausible, and these can also be directly inferred from 
the Balmer decrement with data of sufficient sensitivity. 

5. SEQUENCE OF EVENTS: THE BIRTH AND 
EVOLUTION OF GALACTIC WINDS 

We conclude with a "big picture" model of the star- 
burst wind phenomenon. Our purpose is to provide a 
physical explanation for the delay before the galactic 
wind gets going. We believe that this simple model ex- 
plains the broadest features of our new SPIRAL obser- 
vations. Therefore the stages of evolution for a galactic 
wind are as follows: 

(1) Stars are born in dense molecular clouds whereupon 
the hot massive stars produce a strong UV radiation field 
(Lada & Lada 2003); 

(2) The mean free path of Lyman continuum pho- 
tons is greatly extended by stellar winds (and radiation 
pressure) that blow bubbles into the surrounding ISM. 
(Good discussions of this process are to be found in Cas- 
tor eta^. 1975; Weaver a/. 1977; MacLow & McCray 
1988; Tomisaka & Ikeuchi 1988.) The intense UV radia- 
tion ionizes any pre-existing intracloud medium but this 
may have been cleared out by an earlier wind. 

(3) The over-pressured warm gas evaporated from the 
surfaces of dense clouds helps to fill the wind-blown 
"voids" and intracloud region. The connection between 
the evaporating cloud skins and the intracloud wind ma- 
terial may help the wind to entrain denser gas from these 
same clouds (Coopered a/. 2008). 

(4) The bulk of the supernova ejecta arises after 8-10 
Myr (Fig. [29) Leitherer et al. 1992). As noted by others 
(e.g. Chevalier & Clegg 1985), the ejecta ultimately ther- 
malizes the diffuse material which escapes as a wind; any 
shock heating of the denser medium escapes as infrared 
radiation (Strickland & Stevens 2000). The UV intensity 
declines dramatically since many of the high-mass stars 
have now exploded. 

Thus we are led to an important conclusion. The large- 
scale wind cannot become established until vast quantities 
of warm diffuse gas are dislodged from the surfaces of 
dense clouds by UV evaporation. The ionised material 
with a sufficiently long cooling time is then mechanically 
heated by the supernova mass/energy input. It is there- 
fore not surprising that, by the time the wind gets going, 
most of the short-lived massive stars have faded or dis- 
appeared. We would therefore not expect to see a strong 
signature of stellar photoionisation in the extended wind 
filaments if the starburst phenomenon arises from a single 
burst of star formation. Our new observations may pro- 
vide the strongest evidence yet that starbursts truly are 
impulsive events. The fact that the signature of AGN 
photoionisation is clearly visible in the wind filaments 
argues that the accretion-disk UV timescale is consider- 
ably longer than the starburst timescale. Furthermore, 
unlike the starburst, the mechanical energy input must 
have occurred early in the life of the AGN, presumably 
at a time when jets can also be launched. We argue that 
our observations are a testament to this simple picture. 

Conversely, if there are starburst sources where stellar 
photoionisation does dominate at large radial distances 
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(cf. ShopbeU & Bland-Hawthorn 1998), we anticipate 
that the circumnuclear regions in these sources are char- 
acterised by an extended episode of star formation (i.e. 
continuous star formation) due to an extended period of 
fuelling and accretion (e.g. sustained merger activity). 

6. SUMMARY 

The role of mass/energy feedback and recycling in 
galactic winds over cosmic time is very poorly under- 
stood. In particular, we are far from a detailed under- 
standing of how gas accretion to the circumnuclear re- 
gions fuels local star formation and nuclear activity. So 
what can we say? 

The galaxy wind phenomenon is well established and 
its influence can be detected in the highest redshift 
sources (Heckman a/. 1990; Strickland a/. 2004b; 
VeiWeux et al. 2005). Galactic winds are associated with 
both star formation (e.g. starbursts), nuclear activity 
(e.g. FR I radio sources), and a combination of star for- 
mation and nuclear activity (e.g. Seyfert galaxies). But 
this raises many questions. How do winds get going, 
how energetic are they, how long do they last, and what 
is their role in galaxy evolution? Do we arrive at simi- 
lar answers for AGN-driven winds and starburst-driven 
winds? 

Here we are able to shed light on the first question 
(S|5]). We have carried out the most detailed and ex- 
tensive survey to date of extended plasma gas in active 
galaxies. We constructed ionisation diagnostic diagrams 
to study the nature of the emission in ten galaxies. We 
find that the wholly starburst-powered sources have ex- 
tended emission that is largely consistent with shock ion- 
isation. This sets an important lower limit on the time 
at which the starburst phenomenon began (< 50 Myr). 
At this time, the stellar UV ionisation sources have long 
since faded from view which essentially rules out a con- 
tinuous star formation model for the nucleus, at least in 
our limited subset of starburst galaxies. An interesting 
question is whether we can identify those sources that 
are in an early starburst phase and before the onset of 
the large-scale wind. 

The second question presents difficulties for us, at least 
in the present analysis. In all of our AGN sources, there 
is evidence of circumnuclear star formation: the present 
data are unable to establish which source of activity is 
responsible for the wind. In future studies, it may be pos- 
sible to shed light on this question in one of two ways. 
X-ray spectroscopy can, in principle, establish the [a/Fe] 
abundance ratio of the outfiowing gas. Enhanced a abun- 
dances implicate the presence of SN II, and therefore the 
central importance of the starburst in driving the out- 
flow. Such an effect may hav e already been seen in dwarf 
starburst galaxies ( Martin etal.^ (2002) , Ott, Walter & 
Brinks (2005)), but we are unaware ot similar studies in 
Seyfert outflow sources. This is an important avenue for 
future research. 

An alternative approach is to carefully trace the pres- 
sure changes over the circumnuclear regions down to 
small galactic radii. O'Connell & Mangano (1978) no- 
ticed that the gas densities, as measured from the [SII] 
doublet, appear to decrease along the minor axis of 
M82. Chevalier & Clegg (1985) interpreted this obser- 
vation as a change of wind pressure in a starburst-driven 



wind. This idea has been explored extensively by Heck- 
man efaZ. (1990) who used the [Sii] doublet to measure 
the electron densities as a function of radius in a subset of 
starburst galaxies. In M82, high pressures have b een in- 
ferred from the [S III] doublet (Houck etal. (1984), Smith 
etal. (2006)), from fine structure lines in the far infrared 
(Lugten et al. 1986; Duffy eiaZ. 1987), and from the non- 
thermal radio continuum (Schaaf eia/. 1989). After cer- 
tain assumptions, pressures can be inferred by combining 
observations of radio recombination lines and the asso- 
ciated free-free emission (Rodriguez & Chaisson 1980; 
Seaquist, Bell & Bignell 1985). The intrinsically higher 
resolution of linked radio receivers should allow the pres- 
sure profile to be traced to the smallest VLBI spatial 
scales. Naively, the pressure profile can be relatively 
steep closer to the black hole compared to an antici- 
pated flatter profile across the energy injection region of 
a starburst (Chevalier & Clegg 1985). Thus, we may be 
able to distinguish between AGN-powered vs. starburst- 
powered momentum by measuring the gas pressure as a 
function of galactocentric radius. We leave this to future 
work. 

In the case of the AGN sources, we find that the ex- 
tended line emission is always powered by the accretion 
disk rather than by shock ionisation. Some of these 
sources are clearly associated with central starbursts that 
may well contribute momentum to the outflow in which 
case age-dating the onset of the outflow is more problem- 
atic. In the case of the AGN-driven outflow, these could 
have commenced as recently as 10 Myr or so, although 
if the starburst dominates, the longer timescale is pre- 
ferred. The latter case is interesting because it would 
indicate that the AGN source has a lifetime approaching 
100 Myr, and that the amount of gas converted to energy 
is of order 10^ Mq. 

In the next paper, we combine the ionisation diag- 
nostics with the emission-line kinematics to learn more 
about the energetics of the outflowing gas. The integral 
held spectrograph approach is clearly very powerful for 
studying extended emission line sources. We stress that 
our methods for deriving and comparing ionisation and 
dynamic timescales can be extended to a wider class of 
active galaxies (e.g. quasar hosts), mergers (e.g. ULIRGs) 
and star-forming galaxies (e.g. dwarf irregulars). It will 
be interesting to compare the inferred timescales with 
the star formation histories extracted from circumnuclear 
stellar populations. In high-redshift galaxies where stel- 
lar populations are harder to observe, the ionisation diag- 
nostics may provide key information. In the Appendix, 
we include a brief discussion about how future studies 
can be improved to advance this work further. 
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Fig. 1.— SPIRAL IFU footprints. Images are taken from HST/WFPC2 associations where available and from the DSS for NGC 1482. 
The 22.4" long axis of the SPIRAL IFU is indicated for scale. Note the unusual PA for observations of NGC 6810, NGC 3628 & NGC 5128 
due to the detailed requirements of the NGC 5128 observing program during the May 2007 run. 




Fig. 2. — NGC 1482: Emission line integrated flux maps are genearted by single gaussian fitting to individual spectra. The maps shown 
(left to right, top to bottom) are the continuum at line centre for H« followed by integrated intensity maps (Fint) for Ha, [Olll]A5007, 
[Ol]A6300, [Nll]A6583 & [Sll]A6716. Pixels for which no valid fit was obtained are left blank. 
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Fig. 4. — NGC 1482: The classic lonisation Diagnostic Diagrams (IDDs, see section [2.4| are constructed from the integrated line flux 
measurements. Only IFU mosaic elements with valid line fits for all four (five for [S ll]) Imes are shown in each case. For each set of line 
ratios in turn, classification zones are defined and used to reconstruct the spatial distribution of the corresponding gas in the lower figures. 
The curved solid line marks the extreme starburst limit of Kowlcy etal. (2001). Two fiducial indicators are drawn (straight solid lines). 
The right hand line traces the locus of points for the shock exc ited emission of NGC 1482 (this galaxy). The left indicator traces the locus 
of line ratios for the AGN excited emission of NGC 1365 (Fig. |18|. The dashed line marks the bisector between these two fiducial traces. 
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NGC 253: Emission line maps are generated as for Fig. [2] 




Fig. 6.— 



NGC 253: Line ratio maps are constructed following Fig. [3] 




Fig. 7.— 



NGC 253: lonisation Diagnostic Diagrams are constructed following Fig. [4] and using the same fiducial indicators. 




Fig. 8. — NGC 1808: Emission line maps are generated as for Fig. p] Observations of NGC 1808 were undertaken during a period of 
poor atmospheric transparency and so emission line maps suffer reduced sensitivity. No [O l] emission was recovered. 
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Fig. 9. — NGC 1808: Line ratio maps are constructed following Fig. Is] Observations of NGC 1808 were undertaken during a period of 
poor atmospheric transparency and so emission line maps suffer reduceasensitivity. 




Fig. 10. — NGC 1808: lonisation Diagnostic Diagrams are constructed following Fig.|4] No [Ol] emission was recovred and so the third 
IDD cannot be constructed. 




Fig. 11. — NGC 3628: Emission line intensity maps following Fig. Observations of NGC 3628 were undertaken during a period of 
poor atmospheric transparency and so emission line maps suffer reduced sensitivity. No [O l] emission was recovered. 
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Fig. 12. — NGC 3628: Line ratio maps are constructed following Fig. |3] Due to the reduce sensitivity of the observations for NGC 3628, 
only the [Nil] /Ho map is recovered. 




Fig. 13. — NGC 6810: Emission line intensity maps following Fig. [5] The [O in] map is limited in extent by the sensitivity of the data, 
and no [O l] map was recovered. 
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Fig. 15. — NGC 6810: lonisation Diagnostic Diagrams are constructed following Fig.|4] No [Ol] emission was recovred and so the third 
IDD cannot be constructed for NGC 6810. 
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NGC 1365: Emission line intensity maps following Fig. [2] 
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Fig. 18.— 



NGC 1365: lonisation Diagnostic Diagrams are constructed following Fig.|4] 




Fig. 19.— 



Circinus: Emission line intensity maps following Fig. [2] 




Fig. 21.— 



Circinus: lonisation Diagnostic Diagrams are constructed following Fig. [4] 
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Fig. 22. — NGC 5128 (CcnA): Emission line intensity maps following Fig. [5] 
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Fig. 23. — NGC 5128 (CenA): Line ratio maps are constructed following Fig. [s] Note the data cube was block averaged 3x3 prior to 
spectral fitting in the two [S ll] lines for the construction of the electron density diagnostic due to the weak signal. 




Fig. 24. — NGC 5128 (CcnA): lonisation Diagnostic Diagrams are constructed following Fig. |4] 



38 



Sharp & Bland-Hawthorn 




Fig. 25. — NGC 5128 (CenA): The ionisation cone due the hard AGN spectrum is revealed mo re clearly in the data after spatially 
re-binning the spectra 3x3. The classification scheme from the [Olll]/H/3 vs. [Nll]/Ha IDD (Fig. |24[ l is repeated here. 




Fig. 26.— 



IC 5063: Emission line intensity maps following Fig. [5] 




Fig. 27.— 



IC5063: Line ratio maps are constructed following Fig.js] 




Fig. 28. — IC5063; lonisation Diagnostic Diagrams are constructed following Fig. [4] 
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Fig. 29. — Ionizing luminosity (black) vs. mechanical luminosity (blue) evolution in two different star forming scenarios: (A) instantaneous 
burst; (B) continuous. In both cases, a Salpeter initial mass function is adopted but with two different upper mass cut-offs; 30 Mq (dotted 
line), 100 Mq (solid line). 



APPENDIX 

THE USE OF INTEGRAL FIELD SPECTROGRAPHS - LESSONS LEARNED 

The use of a double-beam integral field spectrograph has been critical to the success of this project. It provides 
both detailed ionisation diagnostics and kinematics over the extent of the outflowing gas. We have demonstrated how 
these IDDs can be used to classify the nature of the wind outflow for each source, and to define spatial zones in which 
different ionisation sources dominate. Furthermore, we have begun the process of spectral and spatial kinematic and 
ionisation deconvolution for each source, in an effort to disentangle the complexity of each system. 

We conclude that blue sensitivity (below 500nm) is a key factor to many of the diagnostic requirements, in particular 
sensitivity to the H/3 and [O ill] spectral region (with SPIRAL we can not comment on the utility of [O ii] emission 
feature for such work). As a clear shock diagnostic, [Oi] is also of great value to classification of the nature of the 
ionisation source, although in all cases presented in this work, sensitivity to the [O ill] emission still presents the 
primary limitation. This limitation is imposed in part by the reduced throughput of the blue arm of the AAOmega 
system ( Sharp e^aL||20(36 ) but will remain an important consideration for more balanced systems. 

Spectral resolution is critical to disentangling the multiple emission components of many of the sources discussed, 
with line splitting in many of the outfiows only marginally resolved in our R^SOOO spectroscopy (~50 km s~^). One 
must of course also consider wavelength coverage. While many of the fainter lines common to ISM studies may be 
beyond the reach of a survey such as ours (for example the temperature sensitive [N ii] A5756 line is typically rather 
faint) we have demonstrated the added value of an extended spectral range in contrast to the limited c overa ge usually 
practical with Fabry-Perot tunable filter systems such as that employed for example by Veilleux et al. ( 2003 ) . A range 
of source redshifts, even for a nearby galaxy sample, further complicates the choice of wavelength range selection. 

The observations and interpretation presented here have been possible primarily due to the high surface brightness 
sensitivity of the AAOmega-SPIRAL system afforded by the large projected spatial pixel size on the sky. The O.Tarcsec 
on a side square elements provide an individual field of view of 0.49arcsec^, which represents factors of ^2, 5.4 & 12 
increase in contrast to 0.5,0.3 or 0.2arcsec systems, more than accounting for the smaller size of the prime mirror of 
the AAT 3.9m when compared to systems available on the current generation of large optical telescopes. This allows 
SPIRAL to move out of the detector readnoise limit, below which re-binning of adjacent apertures will not recover 
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Fig. 30. — Predicted evolution of Ha surface br ight ness of optically thick gas at 1 kpc from a central active source with Lbol = IO^^Lq. 
The black and blue curves are described in Fig. |29| The red curves correspond to 4 different sets of (CA,max;*A) in our AGN model, 
i.e. (0.1, 20), (0.2, 10), (0.4, 5), (0.8, 2.5) where the second value in each pair is the exponential timescale in Myr; larger values therefore 
produce slower decays. The predicted values are at least an order of magnitude higher than what we observe as we discuss in §4. We do not 
include the evolution of the AGN's mechanical luminosity for comparison because this is highly uncertain. But since AGN photoionization 
seems to dominate over shocks, we must assume that the mechanical luminosity profile falls below the AGN UV curve at all times and 
therefore does not resemble the starburst profile. 

signal-to-noise quickly. As recently demonstrated by Westmoquette etal. 09a/b, in truth one must ultimately combine 
the wide field of view and low surface brightness observations practicable with large pixels, with dedicated narrow field 
observations at higher spatial resolution of critical regions of each object. 



